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1. Introduction. 

The topic of finite temperature free field theories on spheres has recently reap- 
peared specifically in connection with the Verlinde-Cardy relation. In a previous pa- 
per, while examining this question, several technical points arose which the present 
work intends to address. The equivalence of the different ways of approaching the 
theory results in, or is a result of, various mathematical identities, most of which 
have been around for some time, many of them classic. It was thought useful to 
present this material in a particular physical context. The opportunity will also be 
taken to present some extensions of our previous programme. 

In [1] two related aspects were considered. One was the symmetry under tem- 
perature inversion and the other the 'ellipticisation' of the total internal energy E. 
It was shown, for example, that at certain temperatures, E could be expressed in 
finite terms involving elliptic integrals. 

In the present paper we discuss, firstly, the technical extension to orbifold fac- 
tors of spheres, S'^/F. This analysis allows us to relate several existing mathematical 
calculations. 

All this detailed work concerns just the internal energy, E, and we turn next 
to our major preoccupation which is a study of the other thermodynamic quantities 
such as the free energy, entropy and specific heat. Roughly speaking, these are 
related differentially to E. 

Our objective, admittedly mathematical, is to determine the behaviour of these 
quantities under modular transformations (which includes temperature inversion) 
and to obtain elliptic expressions for them as far as possible. 

We find for the positive derivatives of the internal energy, such as the spe- 
cific heat, that the results generally parallel those of [1] whereas for the negative 
derivatives (integrals), such as the free energy, there are obstructions to tempera- 
ture inversion symmetry and ellipticisation. The nature of these obstructions makes 
contact with some basic modular concepts, such period polynomials, which are not 
pursued here. 

2. Resume of earlier work. 

Our previous paper [1] was restricted to conformally invariant free field theories 
on the space-times ('Einstein Universes') MxS^ and, apart from a few comments, 
we exclusively discussed odd spheres, S^. We will continue to do so here. 



The approach was based on individual treatment of the terms in the mode 
degeneracy, which is just an odd polynomial in the mode label. The complete sphere 
result was then obtained by simple combination of these individual pieces. From 
this point of view there is nothing special about the sphere. Any odd polynomial 
would do. 

Other treatments, e.g. [2,3], yield more specific forms. We begin here with an 
expression given by Chang and Dowker [4] which is somewhat more general than 
just for the sphere. The situation there discussed was a conformal scalar on an 
orbifold factoring ('triangulation') of the sphere, written S^/F, where F is a regular 
solid symmetry group. 

The total free energy at temperature 1/(3 on MxS^/F was given in the form 

^ oo ^ d ^ 
F = En - - V _e±^o"^/3/2a TT t (I) 

° P^m ^^2smh(d.mP/2a) ^^ 

where the d^ are the degrees associated with the tiling group F and a is the radius 
of the sphere. E^ is the zero temperature, Casimir, value. It can be given in terms 
of Bernoulli polynomials. 

The hemisphere corresponds to (ig = 1; there is only one reflecting plane, and 
all the degrees d^ are equal to 1. The upper sign in (1) refers to Neumann conditions 
on the rim of the hemisphere, and the lower sign to Dirichlet conditions. To obtain 
the full sphere periodic value these two expressions are added, and we will do this 
now for all tilings so as to replace (1) by 

2 oo d .. 

For the time being we will work with just two examples; the hemisphere (giving 
the periodic full sphere) and the quarter-sphere (giving the periodic half-sphere). 
The quarter-sphere is a lune of angle 7r/2 and, in this case, d^ = 2, corresponding 
to the two reflecting planes, also di = 2, d^ = 1 (i = 2, . . . , rf). There are conical 
singularities at the north and south poles. 

For these two examples, more speciflcally, 

j._^ 1 Y^ 1 cosh{m(3/2a) 

2^-1/3 ^^ m smh''{m(3/2a) ' ^ ' 



and 



E, 







1 ^^ 1 coth(m(3/a) 



/3 ^ msmh'^~^{ml3/2a 



Another way of writing (3), which is sometimes useful, is, 

ad ^ 1 
° ^ 2d-^{d-l)d'dd ^. m 



2^-2(rf - 1)13 dp ^^ m2 smh''-\ml3/2a) ' 

As in our previous article, [1], we will work with the total internal energy 

E = aE = ad{(3F/d(3), 

— — 1 ^^ f d — 1 d \ 

° ^ 2^ ^1 Vsinh^-^(m/3/2a) ^ sinh^+^(m/3/2a); ' ^ 

and 

- ^ - ^^ ^ / sech^(m/J/2a) c/ - 1 2(d-l) \ 

° 2^+1 ^^ Vsinh'^+^(m/?/2a) sinh^+^(m/3/2a) sinh^-^(m/3/2a)/ ' ^ ^ 

The full sphere expression (5) involves just the summation 

oo 

Qg{x) = y^ cosech ^mx (7) 

m,=l 

and we now show how to reduce (6) employing also the sum, 

oo 

Rf{x) = y^ sech ■'mx. (8) 

771=1 

Starting from the more general sum 

oo 

Sgj{x) = y^ sech ^mx cosech ^mx (9) 

m,=l 

and writing either 



sech mx = 1 — sinh mx sech m,x (10) 

or 

cosech m,x = cosh m,x cosech m,x — 1 (11) 

we find, for example from (10), the recursion 

/ 



■%.f-]:i'^y{^}% 
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which is easily implemented starting from SgQ = Qg and Sqj = Rf- In particular 

and the energies can be put in the form, if d = 2r + 1, 

E2,+i=Eo + ^(2rg, + (2r + l)g,+i) (12) 

and 

^2,+! = E, + ^ (^(4r - 1) Q, + (2r + 1) Q,^, + Y^{-iy+^-^ Q^ - {-ly R,^ 

(13) 
where the arguments of all the Q's and i?'s are (3/2a. 

In these two cases at least, we see that the problem reduces to the calculation 
of the sums, (7) and (8), which are the subjects of the papers by Kiyek and Schmidt 
[5], Ling [6] and Zucker, [7]. Values of the sum, Qg, for specific g are quite old. 

3. Hyperbolic summations. 

The methods in the papers just referred to, involve elliptic functions, in one way 
or another, which is not surprising in view of the partial fraction representation of 
hyperbolic functions and the double sum form of elliptic functions which originated 
with Eisenstein, [8]. 

Going back to the beginning, Eisenstein, following Euler, discussed the single 
sums (see Hancock [9] p. 32 Ex.5), 

{2g,x)= Y: 7^^7^ = ^''E(-l)'^'^2.,2.cosec2'=7rx. (14) 

n= — oo ^ ' fc=l 

Eisenstein says that the coefficients A2g 2k ^^^ simply related to Bernoulli numbers 
and are given by a recursion derived from continued differentiation. 
Transcribing to the hyperbolic case, 

oo 1 g 

{2g,ix)= V = (-l)g7r^^ VA2g,2fcCOsech^^7ra;. (15) 

n= — oo ^ ' fc=l 
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The explicit form of the recurrence is given by Ling, [6] , and no doubt elsewhere 
after one hundred and sixty years. It is 

^2,+2,2fc = 2g{2l+l) ^^^^ " ^^^^^ " ^^ ^^9,2k-2 + 4/c' A2,,2fc) (16) 

and trivially ^2g,2g = 1- 

Now set X = m^/n and sum again over m referring to the definition (7). One 
finds, 



oo 



m,n= — oo \ ' ' ' k=l 

and we have arrived at an Eisenstein series. 

Equation (17) can be used to determine Qj^ recursively in terms of the left 
hand side which we have shown in [1] is given in elliptic function terms. (This is 
a standard result.) Such is Ling's method of finding the sums Qj^ and is a little 
roundabout. A further instalment is given in [10]. 

The energy (12) can then be determined. Again one has to combine the Eisen- 
stein series values and the calculation is actually not so different from our earlier one. 
The only novelty is that there is now no explicit mention of the mode degeneracies. 
They are, however, lurking in the analysis as the following shows. 

An alternative way of arranging the evaluation of the Q is to rewrite the 
definition (7) as a g-series as done by Zucker, [7]. We re-express his approach 
somewhat to fit our requirements. Resummation produces 



«»-^'E("^rOT5s. 



n=0 

oo 

.2 /„ ix2n ^ r 



22c/ J^ ^2n 



(18) 



-^(in--il,-lf)...n)^ 



^ ^ ' n=0 



where q = e~^, and, not surprisingly, we see here the scalar degeneracies making 
an effective appearance. The full sphere degeneracies emerge on combining the two 
terms in (12) using (18). 

Since they do not occur in our earlier work on the full sphere we now discuss 
the sech sums, (8) which follow from a translation of the argument. 

To begin with, quite easily from (14) and (15), 



(2(7, - - zx) = 2^3 y i --=7r2sV(-l)^+3A2„2fcSech2'=7rx, (19) 

^ y' 2 ^ Z^ (i2x + 2n + lY9 ^^^ ' ^9,2k . v ; 



and, corresponding to (17) 

m,n= — cxD \ ' ' ' fc=l 

(20) 
which contains an 'even-odd' Eisenstein series. Reference to the series in Glaisher 
[11] shows that the Jacobi function zn is involved this time. This has the expansion 

.^zn-u = — tana; — 4 > (— 1)"~^ ^— sin2na:, 

n=l 

where x = 'ku/2K and ^ = 2K/n, which can also be obtained from that for .^zsw 
by adding n/2 to x. Here, as usual, q — e~^ — e~'^^ 1^ ^ and -ftT, K' are given by 

1 

K{k) = [ ^ "^^ =, K' = K{k'), (21) 

where k'^ + k''^ = 1. 

Equivalently the g-series can be constructed, as for cosech. 



n + 2g-l\ g2(n+g) 



«,(,)^2^.^(-ir("+2*-M 



n=o ^ " / -^ y 



-,2n 



(-i)"I^E(("^-(*-i«-")M)"t^,.„. 



(22) 



cf Zucker, [7] , who uses the expansion of nc^ to extract the Lambert series. 

As has been mentioned, our previous analysis [1] expanded the degeneracies in 
powers of n^, n being the mode label, and treated each power separately making use 
of classic elliptic facts. This is the same as Zucker's method and he, more or less, 
reproduces the degeneracy expansion obtained, e.g. by Cahn and Wolf, [12], as well 
as the elliptic expansions to be found in Glaisher, [11]. Zucker's general method of 
finding closed forms for the sums Qg{n) when // corresponds to a singular modulus 
is thus entirely equivalent to our own programme on the sphere, [1]. By the same 
token, the energy on the periodic half-sphere, (6), can also be expressed in finite 
terms using the Sgj, reduced to the Qg and Ri. It does not seem possible to repeat 
this statement for the remaining tilings, (2), at least not obviously. 

The rather negative moral of this calculation is that, in order to evaluate the 
summations that occur for the sphere expressions, (2), in elliptic terms it seems 



necessary to re-introduce the degeneracies in one way or another and to treat the 
individual terms, as in our earher work. The conclusion is that one might as well 
employ the mode-degeneracy expressions in the first place, if these are known, 
particularly since, in addition, some manipulation was required to obtain (12) and 
(13). 

Therefore, although the expressions (2) are particularly compact for all the 
tilings, and provide adequately convergent numerical forms, the earlier piecemeal 
approach has elliptic advantages and so we turn to the degeneracy-mode approach 
for the tilings, knowing that an elliptic formulation is possible for the quarter-sphere 
case at least. We would anticipate that this should be true for all the odd sphere 
'periodic' tilings because the heat-kernel expansion terminates, just as it does on 
the full odd sphere. If this is so, then we could turn the calculation around and 
produce exact forms for the 'new' hyperbolic summations in (2). This will not be 
attempted here. 

4. Mode forms applied to the quarter— sphere. 

In this section we detail the eigenvalue form of the thermo dynamical quantities 
and, as a typical example, consider the specific case of the quarter-sphere, which 
has been treated previously [13] in another connection. 

The Laplacian eigenvalues can be written in the form A = up' ja^ with, [4], 



^■m,n = (a + 2n 4- m) , n, m = 0, 1, 

and the degeneracies are, {d> 1), 

m + (i — 2 
rf-2 

The parameter a = [d ^ 3)/2 for Dirichlet (D) conditions, and a = [d — l)/2 for 
Neumann (N). The D and N expressions are added to get the 'periodic' (P) form. 
For odd spheres we set d = Ir ^ \. For (i > 5 we shift the va label by r + 1 
in the D part and by r — 1 in N, so that the eigenvalues and 'periodic degeneracies, 
9m,u-, read, effectively, 

'^m.n =(2?^ + TO-M), ?i = 0, 1,..., m=l,2, ... 
?Ti-|-r — 2\ f m -\- r 
2r-l y ^ V2r-1 



ym,n 1 n.. i I ' 



We are allowed to extend and adjust the range of m by the vanishing of the bi- 
nomials. For d = 3 attention is needed because the label m runs over 0,1,.... We 
continue with d > 5 and for rf = 3 we only have to add the contribution from the 
m = mode at the end. 

The degeneracies are made more explicit by remarking that g^ is polynomial 
in odd powers of m. Precisely, 

r-l 
k=0 

These are of course related to Stirling polynomials, but we will leave them as they 
are for the time being. 

To proceed with the summation, residue classes mod 2 are introduced so that 
m = 21 + p with < / < oo and p = 0, 1, i.e. m is even (including 0) or odd. The 
two p values are treated separately at first. Then 



cu. 



,^^ = 2{n + l)+p + l = 2N + p + l 



where a further advantageous relabelling has been made to the quadrant coordi- 
nates, 

N = I + n J and u = I — n , 

so that m = N + U + P with A^ = 0, 1, . . . and -A^ < i^ < N . Note that for A^ odd 
z/ only takes odd values and for A^ even z/ only takes even values. 

We now write down the particular spectral quantity in which we are interested; 
the internal energy. 



aE^E^EQ+ 5^ 9^^ 



2lo 

ft 7n .n 



1 — Q ™." 
m,n=U ■' 

where the m label is extended to which is allowed because Qq — Q and we in- 
troduced q = e~'^l^, ^ = {2na)//3. The sums are reorganised to allow for the z/ 
independence of the eigenvalues. Firstly 

oo oo 

m,n=0 P=0,1 l,n=0 

and then 

oo oo N oo TV 

E = E Z+ZE • 

/ n=0 ^^*^ i^=~N N^l i^:=~N 

even even odd odd 



This gives 



^ oo TV \ 2(2Af+p+l) 



N N = l v = -N I ^ 



p=0,l \ iv=o I- 

\ even e^'er^ odd odd 

The summation ranges, A^ even and A^ odd, are rewritten by A^ ^ 2N and A^ 
2N + 1 respectively. The relevant 'degeneracies' are then 

N r-l N 

9even= Y. ^2iV+2.+p = J^ C^ J^ {2N + p + 2u)^''+^ 
v=-N k=0 v=-N 

r-l 27V+P/2 

= 5^C,22'=+1 ^ ^/2fc+l 

fc=0 v'=p/2 

r—1 „2fc 

= E FTY (^2.+2(2iV + 1 + p/2) - B,,^,{p/2)) 

k=0 
N r-l N 

9odd= Yl 92N+2.+p+2 = Y''k Yl m + p + 2 + 2u)^''+^ 

u=-N-l k=0 u=-N-l 

r-l 2Ar+p/2+l 

/2fc+l 



and 



k=0 u'=p/2 

r—1 „2fc 

= E ^ (52fc+2(2iV + 2 + p/2) - B,,^,{p/2)). 



k=o " 






Adding up all contributions, we have 

n2 
+ 

((B2t+2(2iV + 1 + p/2) - B2,+2(p/2)) (4A' + p + 1) ^ " 2(4«+p+i) 



^-^o+EEE^x 



p=0,l Ar=Ofc=0 

2(4Af+p+l) 



2(4Af+p+3) . 



+ (g2fe+2(2iV + 2 + p/2) - B,,^^{p/2)){AN + p + 3) ^^_ ^2(4iv+p+3) ) 
Combining p = and p = 1, 



^-^o + E^x 

A;=0 



E (b2.+2((» + 1)/2)^ - B2>+2(0) ''" : ^^g!;' 



-^2A:+2(l/2; 



271^ 



4n 



'2A:+2V-^/^;^ _ 4^ 
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In order to connect this to elliptic functions we want an explicit expansion of 
the Bernoulli polynomials in powers of the summation index n. To this end we 
apply 



Bj{x + 1/2) = 2^-^B^{2x) - B^{x) 

1 n V / 



k=0 

Noting that -B2A:+i — ^^ ^ = I72, ..., and, combining the various contributions, 
the final answer for the internal energy on the D+N quarter sphere (equivalently 
the 'periodic hemisphere') can be cast into the general form. 



r—L 00 2n 



1=0 n=l ^ 



- „2(2n+l) 



(24) 



n=0 



where the constants /; ^ and J^ are 



-i-i 



, - V "Hj /2 + a + 2jV 2,_, 

r-1 



fc=0 



2(/c + l) 



The zero temperature value, E^q, is a generalised Bernoulli coefficient, [4]. 

This form of the energy, for all dimensions d, is then related to the Fourier 
expansions of zs u and ns u. 

In (i = 5 dimensions the answer reads explicitly. 



2n 1 °o „2n 

3 y 



1 ^-^^ Z?T, 1 ^-^^ 

5 5,0 ^ 24 ^^ 1 - o2" 12 ^-^ 1 

n=l 

1 ^ Q2(2n+1) 

- - > (2?i + 1) ^ ,,, ,, 



O^" 12 ^ — ' i — O^"^ 

n=l n=l 



(25) 



n=0 ~ ^ 



which, equivalently from (13), is to be compared with, 

E, = E,^, + 1(7Q2 + 5Q3 + Qi- Ri)- (26) 

Numerical comparison has also been performed as a check of the calculation. 
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As mentioned, in d = 3 dimensions we have to add the m = contribution to 
the result [24] to get the correct answer. Then one finds, 

_ _ 1 °° ^2n 1 oo 2(2n+l) 

E^^E^r, + -Yn^ ^ , +-V(2n+1) ^ ..o ^,, , (27) 

•3 "^'"2 -^^ 1 — q^^ 2 -^^ 1 — Q2(2n+1) ' ^ ' 

n=l n=0 

with ^3 = 11/480. 

Alternatively, according to (13), 

^3=^3,o + ^(3Qi+3Q2 + ^i) (28) 

and, again, agreement can be shown since, from (18), 
1 1 °° 



^2n 



^2n 
n=l 



For comparison, the energies on the full spheres are 

E 3] full =^2 

— 1 

^ 51 full = ^^(^3-^2) 

in terms of the partial energies, e^, defined in [1] as. 



n=l 

As we have mentioned, the computation via the hyperbolic summations is an 
unnecessary detour only in that, to evaluate elliptically, it seems that the degenera- 
cies have to be reintroduced. 

5. Temperature inversion. 

Apart from the final term, the structure of the energy, (24), is similar to that 
on the full sphere and so, for this part, the conclusions will be the same as in our 
previous paper, [1]. The differences are due to the effect of the polar singularities 
which also account for the final term in E. Because of this term, the recursion 
formula for the e^ employed in (1) now shows that any energy can be expressed, as 
a polynomial in £3 and Cg, plus a multiple of this last term. 

12 



The final term in (24) also affects the bevaviour under temperature inversion, 
^ -^ 1/^. To investigate this, it is rewritten 



q2(2n+l) ^ ^2n ^ 



4n 



/-^^ ^ X _ Q2(2n+1) Z.^ \ — Q2n Z_^ ^ — 0'^'^ 

n=0 ^ n=0 ^ n=0 ^ 

and, while the first term on the right has a well defined behaviour under ^ -^ 1/^, 
described in [1], the second apparently does not, as it corresponds to a different 
temperature. The conclusion is that the inversion properties on the sphere do not 
carry over to its orbifold factors. 



6. The specific heat and modular covariant derivatives. 

In [1] it was shown how the classic elliptic recursions allow the internal energies 
on all odd spheres to be determined in terms of just two parameters. The same 
statement holds for the specific heat, as we demonstrate. From now on we revert 
to the full sphere. 

Before Weierstrass, Eisenstein introduced the double series 

oo ^ 



m-i ,mo = — oo 



{z + rriiUJi + m20J2y 



denoting them by (n, z). If n > 2 they are the higher derivatives of the Weierstrass 
^function. For the time being we adhere to Eisenstein. When n equals 1 or 2, 
the sums are not absolutely convergent and Eisenstein defines them by a limiting 
procedure which depends on the periods, u^, uj^- We will not expound this here 
(see Weil [8]). It is similar to the later procedures of Glaisher and Hurwitz. 
The basic Eisenstein function, E,Jyz\ has the power series expansion 

^n W = - + (-1)- E ( ^ _ 1 J ^t(-^l' -^2) Z''-"- (30) 



in terms of the more common Eisenstein series 



00 



G't(c^i,W2)= XI 7 ;2t' (31) 

■m^,m^ \m^U)^^ +7712(^2) 



by 

Gt{^i,^2) =Gt{uJi,uj2), t>2 
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and 

Gi{^i,^2) =Gi{uJi,uj2) -\ . (32) 



Gi is the quantity that behaves homogeneously under modular transformations of 
the periods, e.g. [14]. 

The Gf are related to the partial energies by, 

et(o = (-irf^^G,(i,z/o, (33) 

for all t. 

By paralleling trigonometric theory, Eisenstein derived a number of basic iden- 
tities for the 6^ including, effectively, the first order differential equation for p, some 
years before Weierstrass. This equation is the origin of the recursion used in [1] and 
which is repeated here, 

-■'«-^^ !L?/)(!:3! g(2"2)--«'----«'- *^^- ('^' 

However, Eisenstein also derived the important relation for the derivative with 
respect to a period, (see [8]). 

27ridS^ 

■^— =^3-^1^2- 35 

From this, one can obtain an equation for a 'partial specific heat', a^, defined by 

..^(^)^.. 

in terms of the more convenient quantities. 

Our notation is that, for example, e = e(^) = e(r) = e(/9), depending on whether 
we wish to use ^, r or /3 as the variable. 

By choosing o^j^ = 1, 0^2 = i/i — r — i/3/27ra and substituting the power series 
into (35) one easily finds 

5 
6 

2t + 3 tzl / 9+ \ (37) 



De,=-2el + ^e2 
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displaying a dependence on Ci which is vital for maintaining the proper inversion 
behaviour under ^ -^ 1/^, which is, 

e,(l/a = (-l)*^e,(0, t>2. (38) 

A formal way of describing the situation is to define a 'modular covariant derivative', 

\ (39) 

so that Ve^ transforms under inversion, ^ — ^ 1/^, like e^_,_]^. Note that, as usual, the 
covariant derivative depends on the behaviour of its operand and we shall extend T) 
to any modular form of weight 2t. (The terminology is that the homogeneous form, 
G^{uJi,uj2), is said to have dimension — 2t while the inhomogeneous, holomorphic 
form, (jt(l, r), has weight 2t.) Although not a modular form, Ci can be said to have 
weight 2. Note that Vci differs from the general structure ofVe^. This definition is 
not that used by Lang, [15], who employs a covariant derivative, 9, in a discussion 
of mod p modularity. 

The general modular form of weight k is an isobaric polynomial in the two in- 
variants, (72 and g^, of elliptic function theory. This follows essentially from the re- 
cursion (34). (Equivalent to ((72, g^) are (e2, €3), (6^2? ^3) ^^^ Ramanujan's (M, A^), 
defined later.) Thus the general form, JF, is written 

a,b:4a+6b=k 

In this setting, the combination in V appears in Ogg, [16], pp. 17, 18. Ogg uses 
the modular transformation behaviour of ei to show that T) increases the weight by 
2 but we could in fact determine this behaviour immediately from (37). Hurwitz, 
[17] §4 employs this device. 

Higher derivatives can be constructed simply since the usual rules apply. For 
example (37) and its derivative read, for all t, 

_ 2t + 3 _^/ 2t \ 

'* "2(2t + 1) '*+! ^ [2k - i/'^'^-^+i 

^-2 (40) 

^2 2t + 3 ^ o v^ /^ 2t \ ^ 

^'^=WTT)^'^+'-^2^^[2k-ir'-'^-'+'^ 

and so on. An equivalent set of equations is given by Ramanujan [18] pp. 165, 166, 
but not using a covariant derivative. 
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Before drawing any consequences from these equations, some relevant elliptic 
facts are interpolated. 

The dependence on Ci means that the complete elliptic function of the second 
kind, E, now makes an appearance via, 

Zucker, [7], from results of Ramanujan, says that E, as well as K, can be expressed 
in terms of algebraic numbers and Gamma functions at singular moduli. 

(41) can be transformed by use of the standard differential relation between K 
and E [e.g. Fricke [19] I,p.46) 

2ek'^^ = E-k'^K 
to which we can also add, 

D^\{^Wf4r.^ ^-f . (42) 

which ellipticises D. In fact, direct manipulation with g-series gives, 

ei = -^Dlog(A;A;'J^3), (43) 

which connects with the alternative expression for e^^, 

ei = -L>log?7, (44) 

in terms of Dedekind's ?7-function, 

r7(r)=Qi/i2jJ(l-Q2j), g = e^-r _ g-v? . (45) 

j 

This includes both the Casimir (zero temperature) contribution and the statistical 
mode sum. 

We have thus more or less arrived at Jacobi's result 

^"=(^(^i-27^i) = ^^'^"^"^ (46) 

another proof of which occurs shortly. 
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Returning to (40), it is apparent from the standard recursions, (34), (40), that 
the ordinary derivative, -De^, can be expressed as a polynomial in e^^, £2 and 63 
(linear in e^^, if t > 1) and the conclusion is that, at temperatures corresponding 
to singular moduli, not only the internal energy but also the specific heat can be 
expressed in finite terms in algebraic numbers and gamma functions. Furthermore, 
the same statement holds for all higher derivatives at singular moduli, as follows 
by iteration of (40) and elimination. Incidentally, it is, of course, possible to derive 
(37) in Weierstrassian vein, and some might prefer this. Such a version is detailed 
by van der Pol [20]. 

Kaneko and Zagier, [21], refer to the algebra generated by e^^, 62 and £3 as the 
algebra of quasi-modular forms. 

For the circle, the three-sphere and the five-sphere we find for the total specific 
heats for conformal scalars, 



a 



(1) 



(^) = H- 



r(3) 



a 



(5) 



(0 = H- 



iO 



2% 
1 

27r 

1 

2n 

T 



.2 5 ■ 




-.«2 + 10 


^3 



(47) 



-^1^3 + 



100 2 2 



120 



As a particular numerical case, at the lemniscate point, ^ = 1, (see [1]), one 



finds 



aW(l) 



1 



■ (rS(l/4) + 247r2r'*(l/4) - 2887r'*) ^ 0.380810377 . 



11527r4 

Other examples we leave to the reader. 

Although equation (37) provides a systematic method, another means of finding 
the derivative in (36) is to apply it to the explicit forms of the partial energies, e^, 
in terms of k and K. This is done by Glaisher [11], §§78,79, who gives the necessary 
rules. One can easily see that the general conclusion will be the same. The structure 
of e^ is a polynomial in /c^ multiplied by a power of K. The ellipticised form of D, 
(42), shows that e^ will appear through the action of d/dk"^ on the K factor. 

Yet another method is to write the Eisenstein function, £^, in terms of the 
^]^-function, and then use the differential (heat) equation that this satisfies, or use 
standard theta-function identities. In modern parlance this is an approach via 
Jacobi forms, which are functions of two variables. We remark that instead of 
the 6'-functions one can use the Jacobi elliptic functions, sn u, etc. The physical 
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significance of tlie u variable could be the radial separation of source and field point 
on the sphere {e.g. Allen et al [22]). 

We now make some further remarks on the recursions and their history. Be- 
cause of the recursions, it is sufficient to give the derivatives of just e^^, 63 and eg, 
and, in order to express these with reasonable factors, it is convenient, and conven- 
tional, to renormalise them. We will use Ramanujan's L, M and A^ notation where 
L = -24ei, M = 240e2 and A^ = -504e3. Then from Eisenstein's equation, (37), 
one rapidly finds the relations often employed by Ramanujan, 

DL= — (L'^-M), or VL=-—M 
12^ '' 12 

DM =-{LM~N), or VM = --N (48) 

o o 

D N ^-{LN - M'^) , or VN^--M^. 

Zi Zi 

An important result now follows easily on noting that from (48) 

L'log(M^-A^2) =L, (49) 

and so, from (44), 

D log(M2 - N'^) = 24D log?7, 

whence, again, Jacobi's relation between the ry-function and the discriminant, A = 

7724 = J— (M^ - Ar2) = -i_ A , (50) 

' 1728 ^ ' (27r)i2 ' ^ > 

the 1728 being the coeflficient of g^ in M"^ — A^^. 

The Weierstrassian way of deriving this relation is somewhat more involved, 
e.g. [23], [24]. Even Weil's reconstruction of Eisenstein's approach mimics this p 
method, [8], p. 33. The modular form proof, [24-27], is quite different and depends 
on the uniqueness of cusp forms of weight 12. 

Equation (49) is equivalent to the statement that A is covariant constant, (put 
2t = 12 in (39)), 

r>A = (D-L)A = 0, (51) 

which highlights the special connecting role played by L. This condition also follows, 
as noted by Tuite, [28] , from the fact that there are no (elliptic) cusp forms of weight 
14 and so works the argument in reverse. 

Equation (51) is a disguised version of the statistical mechanical relation, E = 
— d log Z/d(3 , on the circle. 
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Incidentally, the use of D and T> allows one to give a more general and system- 
atic treatment of Ramanujan's Notebook Entry 45 (p. 352), see [29] part V, p. 484. 
Also, as a side comment, we note that equations (48) and (49) allow derivations of 
the nonlinear differential equations satisfied by M, A^ and A, e.g. [30,20], relatively 
easily. 

The relations (34) and (37) were first obtained by Eisenstein. Later derivations 
have been given by Ramanujan, [18], as mentioned, van der Pol [20], Rankin, [31], 
and Skoruppa [32] . Rankin and Skoruppa use a number theoretic approach, aspects 
of which are reminiscent of Glaisher's work. 

7. Elliptic formulation of the free energy on the circle. 

Equation (44) allows one to obtain an expression for a 'partial free energy', fi, 
defined by, 

- ^ (^f 

i.e. 

A = |- log 77, (52) 

On the circle the total (scaled) free energy is the well known equation, 

^1 = 2A - f log^ = ^ log(r7V0 • (53) 

There is a degeneracy factor of 2, and the log^ term is the zero mode contribution. 
Relation (50) now enables one to give an elliptic interpretation to the partial 
free energy, (52), by 

A = J- log(A/c) , c = (27r)i^ (54) 

Thus the free energy on the circle can be thought of as a function, but not 
an algebraic one, of the elliptic modulus, /c, remembering that ^ = K/K' . See the 
explicit form in (43), (46). 

The choice of constants, in (52), is governed by the desire to satisfy Nernst's 
theorem i.e. to make the entropy vanish at absolute zero, ^ = 0. A partial entropy 
on the circle can be defined by Si = (A — ^i)/C = ~dfi/d^ which is again elliptic. 
It should be said that, when we define these 'partial quantities', we are leaving the 
zero mode contribution in (53) aside as it causes problems with Nernst's theorem 
which have to be addressed separately. 
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8. Higher sphere thermodynamics. 

The question is now whether the higher sphere statistical mechanics can be 
'ellipticised' in the same way as on the circle. The internal energy has aheady been 
treated in [1] but the problem is the free energy, and thence the entropy, which 
demands an effective integration. 

It must be remarked however that Zucker, [33], following Selberg and Chowla, 
uses (46) to derive closed form values for K at singular moduli from computed 
values of ry via the Epstein ^-function. If this were the only method, then the 
elliptic character of the circle free energy would be secondary. However there are 
other ways of finding the singular K which do not involve Kronecker's limit formula. 
We can still therefore maintain the attitude of our earlier work, [1], and regard 
the question of the possible ellipticisation of the free energy on all odd spheres as 
motivation for further investigation. This will be undertaken in a further more 
technical communication. 

8. Summary. 

It has been shown that the temperature inversion properties of the internal 
energy on odd spheres does not carry through to their orbifold factors. The specific 
case of the quarter sphere is considered in detail. During the analysis, hyperbolic 
summations arise that have not so far been encountered. 

It is also proved that the specific heat on any odd (full) sphere can be expressed 
in terms of just the three partial energies, e^^, 63 and eg, the novelty being the 
appearance of e^^, related to the discriminant The derivation makes use of the notion 
of modular covariant derivative with e^ as a 'connection'. We note, incidentally, that 
the discriminant is covariant constant. 

Finally we have demonstrated that the free energy on the circle can be ellip- 
ticised but have left the corresponding statement for all (odd) spheres open and 
subject to later analysis. 
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